AD- A 106  126  ATOMIC  ENERGY  RESEARCH  ESTABLISHMENT  HARWELL  (ENGLAND)  F/G  18/10 

THE  POSSIBLE  USE  OF  CERMET  FUEL  IN  THE  DIDO  AND  PLUTO  HEAVY-WAT— ETC  (U) 
AUG  81  T  D  KENNEDY 


UNCLASSIFIED  AERE-R-10270 


AERE  -  R -10270 


ABSTRACT 

\ 

V, 

J’ International  restrictions  on  the  supply  of  highly  enriched  uranium  have  resulted  in 
the  requirement  to  fuel  research  reactors  with  a  lower-enrichment  uranium  fuel. 

A  study  has  been  made  of  the  feasibility  of  using  low-enrichment  fuels  of  a  new  type 
in  the  DIDO  and  PLUTO  reactors.  This  work  has  been  done  as  a  contribution  to  the  studies 
currently  being  carried  out  internationally  on  the  implications  of  using  lower-enrichment 
fuels  in  heavy-water-moderated  research  reactors. 

The  uranium  content  of  the  U/Al  alloy  at  present  used  cannot  be  increased  sufficiently 
235 

to  maintain  the  requisite  U  i  content  without  undesirable  effects  on  the  physical 
properties  of  the  alloy.  A  different  type  of  fuel  will  therefore  be  required  to  maintain 
the  desired  nuclear  characteristics.  A  possible  solution  to  the  problem  is  the  use  of  a 
cermet  (UjOg/Al)  fuel  material. 

This  paper  describes  a  preliminary  study  of  the  feasibility  of  this  concept  from  the 
heat-transfer  and  safety  viewpoints. 
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Introduction 


The  DIDO/PLUTO  fuel  elements  currently  used  are  of  two  main  types.  Both  consist  of 
a  central  cylindrical  aluminium  thimble  surrounded  by  a  nunber  of  concentric  cylindrical 
fuel  plates  and,  on  the  outside,  an  aluminium  cylinder.  The  layout  of  a  Mk  5/7  fuel 
element  is  shown  in  Fig.  1.  The  inner  (Al)  tube/experiment  thimble  has  a  diameter  of 
5.08  cm,  and  is  used  for  housing  in-core  experiments.  The  other  type  of  element,  known 
as  the  S2  fuel  element,  is  similar,  but  has  more  fuel  tubes  and  a  2.54  cm  central  thimble. 
This  paper  investigates  the  thermal  problems  involved  in  replacing  the  U/Al-alloy 
( 75 %  w/w  U  )  fuel  plates  with  U/Al-cermet  (2095  w/w  U  )  plates,  having  the  same  outside 
dimensions  to  retain  the  same  hydraulic  characteristics.  The  investigation  is  limited 
to  the  Mk  5/7  fuel-elanent  configuration.  Cermet  fuel  plates  have  already  been  used  in 
the  Puerto  Rico  Nuclear  Centre  pool  reactor  (PRNC),  and  assumptions  on  cermet  composition 
are  based  on  information  given  by  Kucera  et  al^'^  for  this  reactor,  which  is  a  5  MW(th) 
pool  reactor.  The  interface  conductance  between  cermet  and  cladding  is  an  unknown 
quantity,  since  this  type  of  fuel  plate  has  not  been  manufactured  in  production  quantities 
in  the  U.K.  Experimental  measurements  will  be  required  if  this  fuel  is  to  be  used;  a 
poor  contact  at  the  interface  could  lead  to  high  temperature  differences  and  possible 
melting  of  the  aluminium  in  the  cermet. 

2.  Fuel-Plate  Dimensions 

From  Ref.  (1),  each  PRNC  fuel  plate  contains  10.67  g  of  U  ,  and  measures  23%  in.  x 

2%  in.  x  0.024  in.  Hence  the  volume  of  cermet  is  1.4175  in3  (23.229  cm3).  Each  DIDO 
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Mk  5/7  fuel  element  contains  205  g  U  .  Thus,  assuming  that  the  DIDO  cermet  fuel  will 
have  the  same  composition  as  the  PRNC  fuel,  the  required  volume  of  cermet  is  easily  found: 

V,  =  23.229  x  =  446.2  cm3. 

The  fuel  loading  in  the  four  plates  of  a  Mk  5/7  element  has  been  optimised  by  Hopper^ 2\ 
and  the  resulting  dimensions  are  given  in  Table  1  below. 

TABLE  1 
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Dimensions  of  Mk  5/7  Fuel  Tubes  and  U  Loading 


Tube 

U23S  (g) 

^  (cm) 

rQ  (cm) 

r  (cm) 

Thimble 

- 

2.540 

2.699 

2.620 

'A' 

43.8 

3.010 

3.178 

3.094 

'B' 

49.3 

3.487 

3.655 

3.571 

'C 

55.8 

3.965 

4.133 

4.049 

’D’ 

56. 1 

4.442 

4.610 

4.526 

Outer 

_ : _ 

4.837 

4.992 

4.915 

The  fliel  loading  in  Table  1  results  in  the  inner  fuel  tube  ’A’  having  the  highest 
fuel  density,  so  calculations  will  be  restricted  to  this  tube  as  the  worst  case.  In  fact, 
the  variation  in  fuel  density  from  ’A’  to  'C'  is  very  small  (x  2'2%) ,  while  tube  ' D ' 
requires  -  less  fueL  density. 
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k  ‘rtmuiv  4*  . 


The  Mk  5/7  fuel  is  'tapered'  at  top  and  bottom, so  that  the  effective  length  is 
18.5  in.  (46.99  cm)  instead  of  24  in.  (60.96  cm). 

Let  the  thickness  of  the  cermet  fuei  in  the  tube  be  t  an.  For  fuel  tube  'A' , 

V  =  23.229  x  =  95.354  cm3. 

C  IO467 


=  2X  rA  It,  so  t  =  Vc/27t  rA  1  = 


95.354 _ 

2 %  x  3.094  x  46.99 


=  0. 1044  cm. 

To  retain  the  existing  'A'-tube  dimensions,  this  would  allow  0.0318  cm  (0.0125  in.) 
for  the  A1  cladding.  For  tubes  ’ B '  and  'C'  the  difference  is  very  small,  while  for 
tube  'o'  the  fuel-core  thickness  required  reduces  to  0.091  cm  (0.036  in.),  allowing 
0.038  cm  (0.015  in.)  for  the  cladding.  Table  2  shows  the  radial  dimensions  of  the  cermet 
fuel. 

TABLE  2 

Dimensions  of  Cermet  in  Fuel  Tubes 


Tube 

r?  (cm) 

r2  (cm! 

'A' 

3.042 

3. 146 

'B' 

3.520 

3.622 

'C' 

3.998 

4. 100 

'O' 

4.480 

4.572 

r(  is  the  inner  radius  and  r 2  the  outer  radius  of  the  cermet  in 

eaclttube 

3.  Steady  Conditions  for  the  Highest-Rated  Fuel  Elanent  at  1.3  x  Nominal  Power  (  1.63  MW) 

Approximately  10$  of  the  fission  heat  from  a  DEDO/PLUTO  MTR  fuel  element  is  deposited 

outside  the  fuel  tubes,  so  the  heat  generation  in  the  fuel  tubes  is  *  1.5  MW  for  a  total 
(2) 

heat  output  of  1.63  MW  .  Fig.  2  shows  the  heat-generation  rate  in  fuel  tube  'A'  of  a 
typical  fuel  element  (03  in  DIDO)  as  a  function  of  axial  position.  The  maximum  power 
density  is  seen  to  be  3020  W/cm3.  The  fuel  tubes  are  sufficiently  thin  that  they  may  be 
treated  as  flat  plates  without  significant  loss  of  accuracy  when  calculating  temperature 
distribution. 

Temperature  difference  from  fuel-centre  plane  to  fuel/cladding  interface 

T  i  .  \  2 


I( t/2) 


where  t  =  fuel-plate  thickness 


cermet  thermal  conductivity  (Fig.  3;  Appendix  1) 
(36.6$  v/v  l^Og) 

0.104  cm,  kc  =  0.32  W/cmK  ; 


AT  = 
c 


3020  x  0.052 
2  x  0.32 


3.1  Temperature  Difference  Across  Fuel/Cladding  Interface  (ATinfc) 

2 

The  interface  conductance  (hint)  has  been  calculated  (Appendix  2)  as  1.42  W/cm  K. 
Ag  =  2*  (Tj  +  r2)  =  2 X  (3.042  +  3.146)  =  38.88  cm2/cm, 

where  Ag  =  interface  surface  area/unit  length. 

Maximum  linear  paver  rating, 


Vx  =  w (*  (r2  -  r?» 

=  3020*  (3.1462  -  3.0422)  =  6  105.8  W/cm 


AT.  . 
int 


Snax 


6  105.8 


A  .h.  „ 
s  int 


38.88  x  1.42 


=  110.6  K. 


3.2  Temperature  Difference  Across  Cladding  (AT  ^) 


AT 


0.t 


cl 


cl 


"cl 


(2) 


(3) 


cl 


=  a  =  =  157.0  w/cm2  =  1.57  MW/m2 

A  38 « 88 


=  0.032  cm 


Kcl  =  2.4  W/cm  K. 


AT 


157  x  0.032 


=  2.1  K. 


Cl  2.4 

3.3  Temperature  Difference  Across  D^O  Coolant  Film  (ATf) 


From  the  subcooled  boiling  correlation  of  Thom 


(3). 


<TW  -  Tsat)  =  22.65  0°* 5  exp(-  P/87)  , 


(4) 


where  T 


w 


wall  temperature  (  C) 


T  =  saturation  temperature  of  fluid  ( 104  C) 
sac 

2 

0  =  surface  heat  flux  (1.57  MW/m  ) 

P  =  fluid  pressure  (1.1  bar). 


T  =  104  +  22.65  /T757.exp(-  0.01264)  =  132.0  C. 

w 

The  D20  temperature  will  be  50°C  at  inlet  and  up  to  74°C  at  outlet,  so  equation  (4) 
gives  a  max.  heat-transfer  coefficient  (h^,)  for  the  D20  at  74°C. 

hf  =  p/ATf  =  1.57/(132-74)  =  0.0271  MW/nA  *  2.71  W/cm2K. 


As  this  coefficient  is  of  the  same  order  as  a  single-phase  coefficient  to  D20  at  high 
mass  flux,  the  latter  will  be  calculated  from  the  Dittus-Boelter  correlation: 

(Nu)  =  0.023  (Re)°,8(Pr)0,4  (5) 
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.  J 

i, 


Total  flow  through  fuel  element,  mT  =  17,300  g/s 

D20  flow  area,  A,,  =  <(13.012  -  2.6992)  +  (3.4872  -  3.1782)  +  (3.96S2  -  3.6552) 
+  (4.4422  -  4.1332)  +  (4,8372  -  4.6 12)] 

=  34.53  cm2 

Mass  flux,  G,j,  =  m^/Ap  =  17,300/34.53  =  501.0  g/an2  s 


(Re)  = 


Grde 


d  = 
e 


4Ar 


4  x  34.53 

2TC(Z(rj+r0)) 


=  0.5782  an 


M- 

k 


=  5.5  x  10"  g/an  s  (viscosity  of  D^O  at  60  C) 


=  0.0063  W/an  K  (thermal  conductivity  of  D20  at  60  C) 


(Re)  = 


501  x  0.5782 


5.5  x  10 


.-3 


=  52,669 


<Pr)  =  3.6;  •••  (Nu)  =  0.023  x  52, 669°‘8  x  3.6°-4  =  229.9 


hf 


(Nu)  .k  229.9  x  0.0063  - - ,  2„ 

=  “d  o75782  =  2‘505  W/cm  K‘ 

e 


•  The  boiling  coefficient  (Thom),  hf  =  2.71  W/cm  K,is  slightly  higher  than  the 
single-phase  coefficient  for  D20  at  the  maximum  temperature  of  74  C,  and  the  maximum 
wall  temperature  of  132°C  is  therefore  selected. 

3.4  Overall  Temperature  Difference,  Fuel  Centre  to  P20  (AT^) 

ATt  =  ATc  +  ATint  +  ATcl  +  ATf  =  12.0  +  110.6  +  2. 1  +  58.0  =  182.7  K 


Max.  fuel-centre  temp.,  T^  =  Tf  +  AT?  =  74  +  182.7 

3.5  Fuel  Mean  Temperature  (at  max.  fuel-temp,  level) 


256 . 7  C. 


t/2 


TU0  =  -i- 


Ix 

TU0  "  2kuo 


•  dx 


t/2 


vo  6\o 


=  256.7  -  (3020  x  0.0522)/(6  x  0.34)  =  252. 7°C. 

Effect  of  25  MJ  Transient  in  Whole  Reactor  Core 

4.1 

In  this  transient,  which  is  required  to  be  withstood  by  the  Harwell  MTRs ,  it  is 
assumed  that  25  \U  of  heat  is  added  to  the  fuel  plates  instantaneously*  ,  i.e.  with  no 
heat  transfer  to  the  coolant.  The  initial  conditions  are  for  normal  full-power 
operation.  Thus,  the  conditions  may  be  taken  from  the  previous  section  for  the 
inner  fuel  tube  'A'  of  a  fuel  element  operating  at  1.63  MW  (  1.3  <  nominal  full  power 
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in  the  highest-rated  fuel  element).  A  25  MJ  whole-core  transient  would  add  1.2  MJ 
to  the  highest-rated  fuel  element  and  (from  the  paver  distribution  with  the  element 
given  by  Hopper^2^)  247.4  kj  to  fuel  tube  'A'. 

Fuel  Tube  'A' 

Volume  of  cermet,  Vc  =  95.354  cm3 

Mass  of  U233  =  43.8  g  (Hopper^) 

Mass  of  U238  =  43.8  x  4  =  175.2  g 

Total  mass  of  U  =  43.8  +  175.2  =  219.0  g 

Total  mass  of  UjOg  =  219  x  (3  x  237.4  +  8  x  16 )/( 3  x  237.4) 

=  258.4  g 

Total  mass  of  UjOg/Al  cermet  (64%  w/w  UjOg)  =  403.7  g 

Density  of  cermet,  pc  =  403.7/95.35  =  4.23  s/at? 

Specific  heat  of  U-Ofl,  C_  =  0.316  J/g  K 

UO 

Mass  of  A1  in  cermet,  m..  =  403.7  -  258.4  =  145.3  g 

A1 

Specific  heat  of  Al,  CL  =  0.985  J/g  K 

Al 

Heat  capacity  of  cermet,  (mCp)c  =  (mCp)l;0  +  (mCp)^  =  224.8  JA 

Volume  of  tube  'A'  ,  V^,  =  Tt(3.l782  -  3.01 2)  x  60.96  =  199.09  cm3 

Volume  of  Al  cladding,  =  Vj  "  vc  =  199.09  -  95.35  =  103.74  cm3 

Density  of  Al  cladding,  =  2.7 1  g/cm3 

Mass  of  Al  cladding,  mci  =  V  ^  =  103.74  x  2.71  =  281.1  g 

Heat  capacity  of  Al  cladding,  (mC  )  ,  =  m  ,C  =  276.9  J/K 

P  CX  Cl  P^ 

Heat  capacity  of  fuel  tube,  (mCp)T  =  (mCp)c  +  (mCp)c^  =  501.7  J/K. 

Assume  that  heat  is  uniformly  distributed  through  the  plate  and  cladding: 

247  4  x  lO3 

Temperature  rise,  AT  =  Q/(mCp)^,  =  - -  = 

However,  10%  of  the  fission  heat  is  deposited  outside  the  fuel  plates,  reducing 
the  temperature  rise  to  0.9  x  493.0  =  443.7  K. 

The  maximum  fuel  temperature  at  the  start  of  the  transient  was  previously 
calculated  as  256. 7°C,  and  the  mean  fuel  temperature  at  the  maximum- temperature  axial 
level  was  calculated  as  252. 7°C.  Thus  the  maximum  fuel  temperature  after  the 
transient  is  calculated  to  be  (256.7  +  443.7)  =  700. 4°C.  However,  this  is  40°C 
above  the  melting  point  of  Al  (660. l°C),  so  the  Al  would  reach  660. 1°C  and  start  to 
melt. 

Latent  heat  of  fusion  of  Al ,  =  401.3  J/g 

Hear  to  melt  Al  in  cermet,  Qfug  =  mA1  \  =  145.3  x  401.3  =  5.83  x  io4  J 
Heat  available  to  melt  Al  (Qavj)  is  found  from  the  apparent  temperature  excess 


above  the  melting  point  multiplied  by  the  heat  capacity  of  the  cermet. 

=  (700.4  -  660.1)  x  224.8  =  9059  J 

Qavl  9059 

Proportion  of  A1  in  cermet  melted  =  — -  =  g83jQ  =  0.155. 

yfus 

Thus,  on  the  assumption  of  an  instantaneous  25  MJ  transient  in  the  core,  with 
all  the  extra  heat  uniformly  distributed  in  the  fuel  tubes  (both  cermet  and  cladding), 
15.5%  of  the  aluminium  in  the  fuel  would  melt.  While  this  limited  melting  of  the 

aluminium  in  the  cermet  would  not  constitute  any  safety  hazard,  the  model  on  which  it 
is  based  is  unrealistic  though  generally  pessimistic. 

4.2 

A  more  realistic  transient  calculation  takes  account  of  heat  transfer  to  the 
coolant  during  a  ramped  power  increase,  and  of  the  distribution  of  heat  within  the 
fuel  tube.  To  investigate  this,  a  two-dimensional  ( r-z  geometry)  model  of  a  fuel 
element  was  set  up,  using  the  finite-difference  transient  heat-transfer  code,  2DT^7\ 

The  model  used  is  shown  in  Fig.  4.  It  consists  of  six  concentric  annular 

tubes,  the  innermost  and  outermost  tubes  being  of  aluminium  and  the  remainder  being 

fuel  tubes  ( Al/cehnet/Al  sandwich)  .  D,,0  coolant  flows  upwards  between  the  tubes. 

Heat  transfer  from  the  outside  boundaries  of  the  model  is  ignored;  it  will  have  a 

negligible  effect,  and  any  error  in  the  model  from  this  cause  must  be  pessimistic. 

Steady  conditions  were  computed  for  the  model  with  1.5  MW  total  heat  generation  from 

the  fuel  tubes  and  17.3  kg/s  total  D„0-coolant  flow.  The  distribution  of  power 
( 2)  i 

among  the  fuel  tubes  was  309.2  kW  from  ’A’,  348.2  kW  from  'B',  406.7  kW  from  'C' , 

and  436.0  kW  from  ' D 1 .  The  flow'  distribution  up  the  coolant  channels  was  taken  from 

(8) 

Lorenz  .  The  axial  distribution  of  power  in  the  fiiel  element  is  taken  from  a 
typical  axial  thermal-neutron  flux  profile  in  DIDO  D3  fuel-element  position ^ 10^  -  as 
shown  in  Fig.  2. 

The  maximum  fuel  temperature  (nominal  steady  state)  in  the  inner  fuel  tube  'A1 
was  computed  to  be  247°C  at  3  cm  above  the  core-centre  plane.  Corresponding 
cladding  and  D20  temperatures  were  129°C  and  65°C  (;n  the  coolant  channel  on  the 
outer  side  of  the  inner  fuel  tube  'A').  The  model  steady-state  temperature  distribu¬ 
tion  was  then  used  for  initial  conditions  for  a  one-second  transient  in  which  the 
power  was  increased  exponentially  so  that  an  additional  1.2  MJ  of  heat  was  generated 
in  the  whole  element,  this  corresponding  to  an  additional  25  MJ  in  the  whole  core 
and  to  247.4  kj  in  the  inner  fuel  tube  'A'.  The  computer  program  accepts  heat- 
generation  rates  in  the  form: 

I(r,z)  =  Io(r,z)enrC  ,  where  IQ(r,z)  is  the  heat-generation  rate  at  position 
(r,z)  at  time  x  =  0. 

Assume  that  the  heat-generation  rates  of  all  regions  increase  at  the  same  rate. 

To  find  the  value  of  m  to  give  the  required  total  heat  production,  over  the  Is 
transient, 


-8- 


where  W  is  the  total  power  produced  by  the  fuel  element.  (I(r,z)  will  vary  with 
time  in  the  same  way  as  W.)  In  this  case,  WQ  =  1.5  MW,  Q^,  =  (1.5  +  1.2)  =  2.7  MJ 
(x  =  1.0  s). 

2.7  =  [em  -  1  ] ;  2.7m  =  1.5  [em  -  l] ;  em  =  1.8  m  +  1; 

m  =  1.045868.  I(r,z)  =  IQ( r,z)  exp( 1 .045868  x) . 

This  transient  has  a  doubling  time  of  0.663  s,  giving  a  power  of  4.27  MW  from  the 

fuel  element  at  the  end  of  the  one  second.  Fig.  5  shows  the  transient  power 

increase  and  the  consequent  response  of  the  highest  fuel  temperature  in  fuel  tube 

'A'.  This  is  seen  to  reach  527°C  at  the  end  of  the  transient,  which  is  133°C  below 

the  melting  point  of  aluminium.  A  doubling  time  of  0.663  s  would  be  caused  by  a 

step  addition  of  0.6$  6k/k.  During  the  one-second  transient  postulated  above,  D00 

mean  temperature  would  increase  by  15°C,  which,  from  the  negative  coefficient- of 

reactivity,  would  give  -0.6$  <5k/k;  thus,  to  sustain  the  transient,  an  additional 

0.6%  5k/k  would  have  to  be  added  over  the  one-second  transient.  The  model  assumed 

2 

single-phase  liquid  heat-transfer  coefficients  to  the  D20  of  2.3  to  2.45  W/cm  K 
throughout  the  transient.  In  fact,  the  coefficients  would  increase  at  first  as 
subcooled  nucleate  boiling  increased  (with  increasing  power) ,  but  might  reduce 
drastically  after  about  0.5  s  as  the  bum-out  heat  flux  was  exceeded.  No  bulk 
boiling  would  occur,  and  the  short  time  and  changing  conditions  would  probably 
exclude  film  bum-out.  This  transient  could  not  occur  with  the  safety  and  control 
systems  operated  with  the  Harwell  MTRs. 

The  maximum  credible  addition  of  reactivity  to  the  Harwell  MTRs  is  taken  to  be 
1.2%  5k/k  at  a  maximum  rate  of  0.75%  6kA>  controlled  by  the  vertical  control  rods 
to  give  the  power  transient  shown  in  Fig.  6  (worst  case).  The  resulting  maximum 
fuel-table-’ A'  temperature  is  also  shown.  The  fuel  temperature  peaks  at  283°C,  a 
rise  of  only  36 °C,  insufficient  to  give  any  safety  problems. 

5.  Heat  Transfer  Within  U„0„  Particles 

— - —  — . - .  3— a—  — 

The  accompanying  sketch  shows  the  maximum-sized  U.0g  particle  which  can  be  accommo¬ 
dated  in  the  fuel  plate.  This  is  assumed  to  be  a  sphere  having  a  diameter  of  0.104  cm 
(i.e.  the  full  width  of  the  fuel  plate. 

Al  Cladding 
Cermet 

Al  Cladding 
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T 


Volume  of  sphere,  VUQ  = 
Density  of  UjOg,  = 
Mass  of  U^Og  sphere, 


4  T  -4  ^ 

|  7i  r  =  5.89  x  10  cm 

8.34  g/cm3 

=  (V.p)UQ  =  4.91  x  10" 


g 


UjOg  contains  84.8 %  w/w  U,  so  mass  of  U  in  sphere,  =  0.848  x  4.91  x  io 


-3 


4.16  x  10  3  g 


U  contains  20%  w/w  U233,  so  mass  of  U233  in  sphere,  in2jg  =  0.20  x  4.16  x  10  3 

=  8.33  x  ,0'4  g 


2^^ 

Mass  of  U  in  fuel  tube  'A'  is  43.8  g 

Power  prodited  in  fuel  tube  'A’  is  309.2  kW 

Power  produced  from  sphere  at  mean  tube-'A'  rating  is 

309.2  x  103  x  8.33  x  iO-4/43.8  =  5.88  W 

Max. /mean  axial  rating  is  1.20 

Power  from  sphere  at  max.  tube-'A'  rating  =  5.88  x  |. 20  =  7.08  W 

Heat-generation  rate  in  sphere,  L,n  =  - — —r~  =  1.20  x  io4  W/cm3 

5.89  x  io 


Temperature  difference  from  centre  to  surface  of  sphere, 

2 

at  Wr  _  (  1.20  x  io4)  x  0.0522 

UO  6k^0  6  x  0.018 


300.4  K. 


The  local  heat  flux  through  the  cladding  to  the  D20  from  the  U30g  sphere  will  be 

higher  than  the  average  flux.  Assuming  half  the  heat  passes  radially  inwards  and  half 

radially  outwards,  the  flux,  -  (7.08/2)/(7t  x  0.0522)  =  416.7  W/cm2. 

max 

Temperature  difference  across  fuel/cladding  interface, 

ATint  =  */hint  =  =  293’3IC 


Temperature  difference  across  cladding 


416.7  x  0.032 


ATcl  =  *  tclAcl  = 


=  5.6  K 


(3), 


Temperature  difference  across  D20  film  (Thom  ), 
(Tw  "  Tsat)  =  22>6S  1)0,5  exP(~  P/87> 


=  22.65  x  /77TC7  exp( -0.01 264)  =  45.7  K. 


T  =  1 04  C ;  T  =  149. 7°C 

t  w 

Max.  U  0  temperature,  T  =  T  +  AT  ,  +  AT.  +  AT 

3  8  0°  w  cl  mt  UO 


=  149.7  +  5.6  +  293.5  +  300.4  =  749. 2°C  . 


-10- 


Such  a  particle  size  (1040  nm)  is  much  larger  than  the  upper  limit  (149  pm)  of  the 
range  of  particle  sizes  reported  by  Kucera  et  al^.  Although  this  large  size  does  not 
lead  to  unacceptable  temperatures,  a  more  likely  maximum  size  is  200  pm  (s  0.02  cm). 


This  leads  to: 


1.20  ^  iQ4  x  O.Qi2 

6  x  0.018 


=  1 1.  I  K. 


Such  a  particle  would  not  significantly  perturb  the  average  heat  flux  from  the  cermet, 
so  the  maximum  UT0Q  temperature  may  be  calculated  from  the  maximum  cermet  temperature,  T  , 

O  o  C 

computed  (Section  3). 


Tc  +  ATU0 


=  247+11  =  258  C. 


To  assess  the  effect  of  the  maximum  transient  power  excursion,  illustrated  in  Fig.  6, 
on  a  U-0Q  particle  of  the  size  of  a  200  pm  diameter  sphere,  a  simple  model  was  set  up, 

O  o 

consisting  of  a  right  cylinder  of  UjOg  (174.8  pm  diameter)  surrounded  by  a  hollow  right 
cylinder  of  aluminium  (244  pm  diameter).  The  model  therefore  contained  the  same 
proportion  of  U,0o  (36.6$  v/v)  to  aluminium  (63.4%)  as  the  cermet.  The  U,0o  cylinder 

w  o  0  3 

was  of  the  same  volume  (4.19  x  10  cm  )  as  a  200  pm  diameter  sphere.  The  model  was  set 
up  with  an  all-round  boundary  temperature  of  247°C,  representing  the  general  cermet 
temperature  computed  in  Section  3.  The  U^Og  was  assigned  a  heat-generation  rate  of 
12,000  W/cm3,  corresponding  to  a  fuel-element  power  of  1.5  MW. 


122  pm 


i'.  87-4 pm,  y) 

/  . 

I  ✓  /  /  /  /I 


/  /  / 


\/.t- 


U3  °8 


Computer 
Model  of 
U3  08  Particle 

in  Cermet 


This  model  was  then  run  to  steady  state,  which  showed  a  maximum  U, 0Q  temperature  of 

O  O 

258°C.  From  this  initial  condition,  a  transient  of  the  form  shown  in  Fig.  6  was  applied 
to  the  model;  the  temperature  history  at  the  centre  of  the  UT0„  particle  is  also  shown 

•J  O 


+ti-  ^.4  i,  -  --  r-  •  f— - fr  •/«***>  , 


in  Fig.  6.  The  resulting  maximum  UjOg  temperature  was  297°C  at  0.37  s  after  the  start 
of  the  transient.  The  exponential  transient  I  =  IQ  exp(  1.045868  x))  was  also  applied 
to  the  model  for  one  second,  with  the  results  shown  in  5.  The  maximum  UjOg 
temperature  was  then  582°C  at  the  end  of  the  one-second  transient. 

6.  Conclusions 

Provided  that  good  bonding  can  be  obtained  between  the  cermet  and  the  aluminium 

cladding,  there  appear  to  be  no  problems  from  a  heat-transfer  viewpoint  to  the  use  of 

U^Og/Ai -cermet  fuel.  Where  information  has  been  lacking,  pessimistic  assumptions  have 

been  made  in  the  calculations,  which  nevertheless  lead  to  maximum  cermet  temperatures 

below  any  which  should  cause  problems.  Experimental  measurements  of  thermal  conductivity 

should  be  made  when  the  cermet  is  selected.  Local  temperatures  in  the  U„0Q  particles 

o  o 

within  the  cermet  have  been  predicted,  and  present  no  problems  in  normal  operation,  nor 
for  the  maximum  credible  transient  power  excursion  for  the  Harwell  MTRs. 
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APPENDIX  1 

Thermal  Conductivity  of  U.0Q-A1  Cermet 

Fig.  3  from  Ref.  (4)  shows  the  thermal  conductivities  of  L'jOg-Al  dispersions  at 
ambient  temperature  (20°C). 

The  cermet  considered  has  64$  w/w  UjOg  and  36$  A1 .  Density  of  U^Og  is  8.34  g/cm3 
and  that  of  A1  2.71  g/cm3. 


PA1*VA1  36 


2.71  V. 


.V,„  64  8.34  V,, 


f\jO,VUO 


'  VA1  =  '-73'  VU0* 

The  cermet  has  63.4$  v/v  A1  and  36.6$  UjOg. 

Touloukian^  gives  thermal  conductivities  for  a  number  of  U02-metal  cermets;  three 
of  these  are  shown  in  Table  Al.1  below,  together  vith  the  conductivities  of  UC>2  and  the 
component  metals.  These  cermets  contain  higher  proportions  of  uranium  oxide  than  the 
UjOg-Al  cermet  considered  here.  However,  in  each  case  the  thermal  conductivity  increases 
with  rising  temperature,  although  the  UC>2  conductivity  on  its  own  reduces,  as  do  the 
conductivities  of  two  of  the  component  metals.  U30g  conductivity  is  low  at  ambient 
temperature,  and  would  not  be  expected  to  vary  significantly  with  temperature.  A1  con¬ 
ductivity  falls  fairly  slowly  from  120°C  to  650°C.  On  this  evidence,  it  seems  likely  that 
the  conductivity  of  the  Al-UjOg  cermet  would  increase  with  temperature,  so  an  assumption  of 
conductivity  invariant  with  temperature  is  probably  pessimistic. 

TABLE  Al.1 


Material 

_^2 _ 

Mb 

U02-Mb  cermet 
(70$  v/v  U02) 

St. Steel 


k  j  ( W/cm  K) 

0.0774 

1.371 

0.201 

(0.183)* 


0.079 

(0.094)* 


T,  (°C)  k2  (W/cm  K) 


t2  (°C) 


0.0304 

1.147 

0.204 
(0. 104)* 


0. 142 
(0.058)* 


>  *?  v: 

>  u 

it  oC 


0.980 

U02“Cr  cermet 
(80%  v/v  U02) 

0.146 

(0.129)* 

U3°8 

mmm 

A1 

2.36 

O3O8-AI  cermet 
( 37$  v/v  U30g) 

0.34 
(0. 388)* 

0*018? 


2.13 


r 


A  simple  method  of  estimating  cermet  conductivity  from  its  component  materials  is  to 
take  a  logarithmic  mean  value: 


In  kc  =  a  in  ka  +  3  In  1^  ,  (Al.l) 

where  subscripts  a,  b  and  c  refer  to  the  metal,  the  oxide  and  the  cermet  respectively  and 
a  and  3  ate  the  volume  fractions  of  metal  and  oxide  respectively;  e.g.  for  the  UC^-tMb 
cermets  at  50°C 

In  kc  =  0.3  In  (1.371)  +  0.7  In  (0.0774)  =  -1.6965 

k  =  0.183  W/cm  K  . 

c 

This  value  is  close  (-9%)  to  the  value  of  0.201  W/cm  K  from  Ref.  (5).  Using  this  method 
for  the  U02~Mb  cermet  at  650°C  gives  kc  =  0.104  W/cm  K.  Thus,  while  the  method  provides 
a  good  estimate  of  cermet  thermal  conductivity  at  low  temperature,  it  underestimates  the 
measured  value  by  w  50$  at  the  higher  temperatures.  A  similar,  but  even  more  pronounced 
effect  occurs  with  the  other  two  cermets.  The  values  calculated  from  equation  (Al.l) 
are  given  in  brackets  below  the  values  from  Ref.  (5). 


w 
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APPENDIX  2 


Interface  Conductance,  U„0„-A1  Cermet  to  A1 

“  -it  O 

The  conditions  at  this  interface  are  not  known,  but  because  of  the  very  high  heat 
flux  through  it,  contact  must  be  good  to  prevent  melting  of  the  aiuminium  in  the  cermet. 

(9) 

Main  gives  methods  of  calculating  interface  conductance.  The  starting  point  is  to 
assume  three  separate  conductance  paths  in  parallel: 

hT  =  hf  +  hg  +  hr  ,  (A2.1) 

where  subscripts  T,  f,  s  and  r  refer  to  total,  fluid,  solid  and  radiation  respectively; 
h  is  the  conductance  (or  heat-transfer  coefficient  in  the  case  of  h^,).  In  this  case, 
h  may  be  ignored,  since  any  temperature  above  600°C  at  the  interface  will  result  in 
melting  of  the  aluminium. 

(a)  Fluid  Conductance,  hf 

^f 

h{.  =  -j-  ,  (A2.2) 

where  kf  is  fluid  conductivity,  k^.  (air)  =  0.0004  w/cm  K,  and  d  is  the  effective  gap 
thickness. 


d  =  C(RJ  +  R2)  +  g,  +  g2  , 


(A2.3) 


where  R(,  R2  are  surface  roughnesses,  C  is  the  roughness-height  factor  (C  =  1.5  for 
heavy  contact),  and  g(  and  g2  are  temperature-jump  distances. 


Here  assume  R(  =  R,,  =  25  x  10  cm. 


From  Ref.  (9) , 


8,  -  82 


(2  -  a)Pa 

a  P 


Z  , 


(A2.4) 


and  Z  s  3  x  10~  cm,  a  x  0.45,  P 

-5 


P„  =  1  bar. 

a 


8,  -  82 


1.55  x  3  x  iQ 
0.45 


=  1.03  x  10  cm 


d  =  1.5  (2.5  x  10"6  +  2.5  x  )0-5)  +  2  ( 1 .03  x  10-4) 


=  2.81  x  10  am 

k„ 


hf  "  d  " 


-4 

i 

4  x  10” 


=  1.42  W/cm  K  . 


2.81  x  10 
(b)  Solid  Conductance  h 


This  is  very  difficult  to  quantify,  and  since  the  value  of  by  found  above  is 
adequate  to  prevent  an  excessive  AT  across  the  interface  no  attempt  will  be  made  to 
do  so.  For  the  U-Al-alloy  fuel,  the  assumption  is  normally  made  that  contact  is 
complete,  and  hence  hg  =  «>;  in  this  case,  the  pessimistic  assumption  of  no  solid 
contact  is  made,  so  that  hg  =  0. 
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Nomenclature 


A 


s 


\ 

C 

CP 

d 

de 

g 

G 

h 

X 

k 

1 


m 

P 

Q 

r 

R 

S 

t 

T 

V 

W 

z 


Greek 

a 

AT 

X 

4 

i 

x 


Surface  area 

Total  flow  area  through  fuel  element 

Rougiiness  height  factor 

Specific  heat 

Effective  gap  thickness 

Equivalent  diameter  of  coolant  channels 

Temperature- jump  distance 

Mass  flux  of  coolant. 

Heat-transfer  coefficient 

Heat-generation  rate 

Thermal  conductivity 

Fuelled  length 

Mass 

Pressure 

Heat 

Radius 

Rougiiness  height  of  surface 

Wetted  perimeter 

Thickness 

Temperature 

Volume 

Power  from  fuel  element 
Temperature- jump  parameter 


Accommodation  coefficient 
Temperature  difference 
Latent  heat  of  fusion 
Dynamic  viscosity 
Heat  flux 
Time 


Dimensionless  Groups 
(Nu)  -  Nusselt  No.  (=  hdg/k) 

(Pr)  -  Prandtl  No.  (=Cp.|Vk) 

(Re)  -  Reynolds  No.  (=G.de/px) 

Subscripts 


Ai 

-  Aluminium 

int 

-  Interface 

avl 

-  Total  available 

max 

-  Maximum 

c 

-  Cermet 

0 

-  Outer 

cl 

-  Cladding 

5 

-  Surface 

e 

-  Equivalent 

sat 

-  Saturation 

f 

-  Fluid  or  film 

T 

-  Total 

fus 

-  Available  for  fusion 

UO 

'  ll03°8 

i 

-  Inner 

W 

-  Wall 
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core -centre  plane 


Thermal  conductivity  (W/cm  K) 


T 


0.5 


Uranium  loading  (g/ctc  ) 

1.0  1.5  2.0 


2.5 


3.0 


i  i  i  i  i  i  r 


2.5 


2.0 


1.5 


1.0 


0.5 


Q-Mc Elroy  (1970) 


Q-U’illiams  and  Graves  (1980) 


0.1  0.2  0.3  0.4 

U_0o  in  Core,  Volume  Fraction 

■S  O 


Fig.  3  Thermal  Conductivity  of  U,0_-Al  Dispersions 

3  O 


0.5 


L 


HEAT  -  GENERATION  RATES 
IN  FUEL  TUBES  (  kW ) 
STEADY  CONDITIONS 


ZONE  'A' 

•B' 

•r 

■D' 

1  30  7 

346 

40-3 

43-5 

2  52  9 

59  -  6 

69  6 

74  5 

3  55  8 

630 

73-4 

78  7 

4  55  6 

62  5 

73  0 

78  2 

5  50-1 

56-4 

65  0 

70  5 

6  41-2 

46-4 

54-2 

58  6 

DEO 

24-4 

28-4 

30-4 

TOTAL  1307-9  i 

346-9 

404  7 

434-4 

FISSILE  HEAT  -  1493-9  kW 
GAMMA  HEAT  { Al  )  -  6-1  kW 
TOTAL  HEAT  -  1500  kW 


COOLANT  CONDITIONS 
(STEADY  STATE  ,  1500  kW  ) 


No. 

FLOW 
(kg/s ) 

H.T.Coeff 
(  W/ottK 

tin<°c) 

1  . 

2-768 

2  30 

50 

63-4 

2 

3-244 

2-34 

50 

74-0 

3 

3  -719 

2  38 

50 

74  2 

4 

4-152 

2-42 

50 

74  1 

5 

3-417 

2  -45 

50 

65  9 

TEMPERATURE  RESPONSE  TO  CORE  -  POWER  TRANSIENT  (+0-75%  6k  /  k 

REACTOR  TRIP  AT  1-65MW  ) 


